Hamiltonian Flow in Coulomb Gauge Yang-Mills Theory 
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I. INTRODUCTION 

QCD is experimentally well-tested as the theory of 
strong interactions in the perturbative regime where 
asymptotic freedom holds. In the non-pcrturbative 
regime we have collected much theoretical evidence rang- 
ing from a linearly rising potential for heavy quarks to 
hadronic observables which match at least qualitatively 
the experimental values. These non-perturbative results 
were either obtained by lattice methods, QCD-model 
computations and in recent years, to an increasing extent, 
by first principle continuum QCD computations based 
on functional methods. Despite the impressive success 
of these combined efforts there still remains much to be 
understood, both qualitatively and quantitatively. Open 
physics questions range from the mechanism of confine- 
ment and its relation to spontaneous chiral symmetry 
breaking to the properties of QCD at finite temperature 
and density. Progress in these directions can only be ob- 
tained with a combination of the methods at hand, in 
particular, when it comes to the understanding of the 
underlying physics mechanisms. 

This has triggered an increasing interest in non- 
pcrturbative studies of continuum Yang-Mills theory and 
full QCD in recent years. Many of these studies were car- 
ried out in Landau gauge, using Dyson- Schwinger equa- 
tions (DSE), for reviews see (see also @), as well as 
the functional renormalization group (FRG) equations, 
for reviews see (see also Another attractive 

possibility is Coulomb g a.uge, w hich has been pursued 
either with DSEs, e.g. |ll| - [l3| . or with a variational 
approach to the Hamiltonian formulation p^ - fioj . for a 
short introduction see [2^. Each approach has its own 
advantages and drawbacks and by combining the dif- 
ferent approaches one can expect to gain new insights 
into the theory, in particular, into the non-perturbative 
regime. 

In this spirit we put forward in the present paper a 
functional renormalization group approach to the Hamil- 



tonian formulation of Yang-Mills theory. A specific ad- 
vantage of the Hamiltonian formulation is its close con- 
nection to physics as demonstrated in various applica- 
tions, see e.g. [2l|. However, Hamiltonian renormaliza- 
tion is a well-known non-trivial task which has hampered 
progress in the Hamiltonian approach for many years. 
Specifically it complicates the search for RG-invariant, 
self-consistent approximations to the full vacuum wave 
functional. In the present Hamiltonian FRG approach 
to Coulomb gauge Yang-Mills theory the latter task is 
directly solved by a self-consistent approximation to the 
flow equation. Such a procedure guarantees RG invari- 
ance by its very definition, and hence combines the ad- 
vantages of the Hamiltonian approach with that of the 
FRG. We will also show explicitly that with suitable trun- 
cations the integrated flow equations become precisely 
the DSEs of the variational approach [l^ with a Gaus- 
sian ansatz for the vacuum wave functional. In our ex- 
plicit computations we will focus on the infrared sector of 
the theory within the ghost-dominance scenario. There- 
fore we only include full momentum-dependent gluon 
and ghost propagators and the bare ghost-gluon ver- 
tex. The latter approximation is justified by the non- 
renormalization of the full ghost-gluon vertex. 

The paper is organized as follows: In section II we 
present the basics of the functional renormalization group 
(FRG) flow equation approach. In particular, we derive 
the flow equation for the Hamiltonian approach to Yang- 
Mills theory in Coulomb gauge. The gauge fixing in the 
scalar product of the wave functional is accomplished by 
the Faddcev-Popov method. In section HI the FRG flow 
equations for the gluon and ghost propagators arc de- 
rived assuming a bare, non-running, ghost-gluon vertex. 
We also extend the uniqueness proof for the infrared scal- 
ing solution in Landau gauge to Coulomb gauge Yang- 
Mills theory. In sect. IV we numerically solve the flow 
equations for ghost and gluon propagators within two 
different approximations which are compared with the 
DSE results. Finally our conclusions are given in sect. 
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V. Some mathematical details of the derivation of the 
flow equations of the Hamiltonian approach as well as 
some details of the numerical procedure are deferred to 
the appendices. 



II. DERIVATION OF THE HAMILTONIAN 
FLOWS 

A. Flow equation for the effective action 

Below we briefly summarize the essential ingredients 
of the FRG approach. The starting point is the finite 
renormalized generating functional of the Green's func- 
tions, 



(1) 



where the subscript "ren" of the measure indicates an 
appropriate renormalization procedure that renders the 
functional integral in ([1]) finite, for more details see Q. 
Here ip and j denote collectively all fields involved and the 
corresponding sources. Furthermore, the scalar product 
j ■ tp includes summation over all indices and integration 
over space-time. Finally, the theory is specified by the 
classical action 5 [1^9] and the renormalized functional in- 
tegral measure. In the functional renormalization group 
(FRG) approach the generating functional is IR regu- 
larized by adding a regulator term ASk to the classical 
action 5*. It depends on a cut-off momentum k and is 
chosen to be quadratic in the fields, 

ASM = \p-Rk-p^\j ^p{^p)Rk{pMp) . (2) 

The regulator function Rkip) is an effective momentum 
dependent mass and has the properties 



lim Rk{p) 
lim Rkip) 

fe/p— i-O 



> 







(3) 



The first condition ensures that Rk{p) is indeed an in- 
frared regulator and suppresses the propagation of the 
modes with momentum p ^ k. The second condition 
implies that the momentum modes with p ^ k are un- 
affected by the regulator and that the full finite renor- 
malized generating functional of the theory at hand is 
recovered as the cut-off scale k is pushed to zero. 

The basic idea of the FRG flow equation approach to 
Yang-Mills theory is to start at a large cut-off scale fc, 
where the theory is under control due to asymptotic free- 
dom, and then let the cut-off k flow to the small momen- 
tum regime, which is non-perturbative. The evolution 
of the Green's functions with the cut-off scale k is de- 
scribed by a flow equation which is obtained by taking 
the derivative of the regulated generating functional 



(4) 



w.r.t. the momentum scale k. We also remark that the 
finiteness of Zk[j] follows from that of Z[i] with the rep- 
resentation 



Zkm=e-^^-^"'^^Z\j] 



(5) 



and the differentiability of Z[j] w.r.t. j. With ([5]) the 
flow of Zk is derived as 



dtZk[j] 



\ 5 ■ (5 , 



where the dot on R stands for the derivative w.r.t. the 
dimcnsionless variable t = Infe/feg. Here fcg is an arbi- 
trary reference scale. The functional derivative 8/5j in 
momentum space is to be understood as 



(7) 



Expressing ([5]) in terms of the generating functional of 
the connected Green's functions, defined in (U), 

we get 



dtWk[j\ = 
where 

S^Wk 



l5Wk ■ 5Wk 1^ • S'^Wk 



Tri?fc 



d<^p ■ 



(27r)° 



Rkip) (27r 



S^Wk 



^ii-p)5jip) 



(8) 



(9) 



By taking derivatives of (|8]) w.r.t. j one obtains the flow 
equations for the connected Green's functions. In prac- 
tice, it is usually more convenient to perform first a Leg- 
endre transform from the sources j to the classical field 



5Wk[]] 



resulting in the effective action 

rk[<f>] = i-Wk[j]+j 



(10) 



(11) 



where iki(f) is given by solving ([TU]) for j. Hence, pT|) 
also implies that 
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5iTk + ASk 



(12) 



By taking the derivative of the effective action w.r.t. 
t = Infc/fco and using its definition (|TT|) and ([U one ar- 
rives at 



(13) 



The second derivative of Wk w.r.t. the currents is the 
connected two-point function, the propagator of the the- 
ory. It is related to the inverse of the second derivative 
of the effective action, 



S^Tk 



(14) 
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Eq. (Hi]) follows directly from and 



■J 



S_ 



(15) 



Inserting (|14l) into ([T5|) we obtain Wetterich's flow equa- 
tion for the effective action [23, 



where the dot denotes the derivative w.r.t. t and 



r 



(n) 



(16) 



(17) 



are the one-particle irreducible n-point functions (proper 
vertices). We have also introduced a condensed notation 
where n stands for the space-time variable, = (j){x„). 
The generic structure of this equation is independent of 
the details of the underlying theory, i.e., of the explicit 
form of the action S[(j)\, but is a mere consequence of the 
form of the regulator term ([2]), i.e., that it is quadratic in 
the fields. By taking functional derivatives of Eq. 
w.r.t. the fields one obtains the flow equations for the (in- 
verse) propagators or vertices. Taking the second func- 
tional derivative we find the flow equation for the two- 
point function, 



^(2) 
fc,12 



^ -Tri?fe 
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.(4) 
k,12 



.(3) f-r{2)t 
k,l 



r^^'[cb] + Rk] ri^^ + (io2) 



r 



(18) 



In psp we have suppressed all cyclic indices, which are 
summed (integrated) over in the trace. The flow equation 
psp for the propagator is diagrammatically illustrated in 

Fig. m 



B. Hamiltonian flows 

In this section we derive the flow equation for the ef- 
fective action in the Hamiltonian approach. Some de- 
tails can be found in Appendix [X] In the Hamiltonian 
approach the generating functional of the static (equal- 
time) Green's functions is defined by 

Z[j] = exp(j • (^)IV^) = J Pt^lV^MP exp(j ■ ^) , (19) 

where {(p\tp) = is the vacuum wave functional. Com- 
parison with Eq. ([T]) shows that in the Hamiltonian ap- 
proach the "action" is defined by the wave functional 



exp(-5M)^|^[<^]P 



(20) 



With this identification Eq. ([T9|) has precisely the stan- 
dard form of the generating functional ([1]), except that 
the functional integral extends over time-independent 
fields, and we can repeat, step by step, the formal ma- 
nipulations of the last section to derive the corresponding 
FRG flow equation, which has the same structure as Eq. 

m- 

In the present work we are specifically interested in 
the Hamiltonian flow of Yang-Mills theory in Coulomb 
gauge. In this case (p stands for the transverse spatial 
components of the gauge field A, 



tijAj 



with tij (p) 



r^2 



Note that in the present Hamiltonian context, we denote 
the contravariant components of all 3- vectors with sub- 
scripts. The wave functional ip{A) is the true Yang-Mills 
vacuum wave functional restricted to transverse fields. 
Implementing the Coulomb gauge in the standard fashion 
by the Faddecv- Popov method, the generating functional 
([TOl) reads 



Z[J] 



VAdct{~dD)\'iP[A]\^ c^j>{J ■ A) 



(22) 



where dct{—dD) denotes the Faddcev-Popov determi- 
nant and 



jjab ^ gabg _^ gfacb^c ^ 



(23) 



is the covariant derivative in the adjoint representation. 
The short hand notation J ■ A also includes the internal 
indices. 



J-A = 



(2^ 



jr(-p)^^(p). 



(24) 



The Faddeev-Popov determinant det(— 9D) is repre- 
sented by an integral over ghost fields c''(x), c°(x). Af- 
ter introducing regulator terms for ghosts and gluons via 
AS'fc, see Eq. ©, the regulated generating functional 
(m of the Hamiltonian approach to Yang-Mills theory in 
Coulomb gauge reads 

= expWfe[J,a,a]. (25) 
The "classical action" S in is given by 

(21) S^-\n\^[A]\^ + j dPx(f{^){-dDY''c\^) (26) 



kdk 



FIG. 1: Flow equation for the inverse propagator (2-point function) T\ 



Here and in all the diagrammatic equations 



below, lines with full circles represent ^F^'^' + R^^ on the r.h.s. of the equation, but only F^^' on its l.h.s. Open circles 

in) 

represent full proper vertices F^ , while the square with a cross represents the regulator insertion Rk. 



and the regulator term is chosen as 

A5fc [A, c,c\ ^ ^A- Rjif, ■ A + c- R^j^ ■ c . 



(27) 



Here we have used again the short hand notation for 
scalar products, e.g. 



A • i?^., • A = 



^Al{~vXM^\^V). (28) 



In pS]) the transversality of the gauge field, ^"(p) = 
ty (p)^"(p), is implied. Put differently, we could have 
multiplied the regulator with the transverse projector. 
Due to global color symmetry and spatial rotational sym- 
metry, the regulators take the form 



(29) 



with the notation p = |p| that we shall use from now on. 
Both regulators are chosen to depend on the same dimen- 
sionless shape function rk{p). Accounting for dimensions 
we put 



RcMp) =9p'^ rk{p) = gRc,k{p) ■ 



(30) 



(see our remarks following Eq. ([50)) concerning the factor 
g included in the definition of Rc,k{p))- In the numerical 
solution we have chosen the rk{p) as given in (|57|) . 

From the regularized generating functional (|25p one de- 
rives the flow equation for the effective action as outlined 
in the previous section for a single field. In the present 
case, for compactness of the notation, it is convenient to 
combine gluon and ghost fields into a superfield 



9? = (A,c,c) , ip = (A, -c,c) . 
Accordingly, we introduce the supersources 

/ ~ ( J, cr, (t) , I ~ (J, — (T, a) 
and supermatrices 



(31) 



(32) 



7^fe = Aia.g{RA,k-,Rc.k-,R, 



'c,k) 



M 



diag(l,-l,-l). 



(33) 



where M figures as metric in the superspace and enters 
the definition of the supertrace, 

STr(. . . ) = Tr(Af . . . ) . (34) 

With this notation the effective action PT|) is given by 

1 



-Wk[Ik 



h 



■ TZkM ■ 



where Wfe[/fe] is defined by Eq. (|25p and 
7 SWkih] 



(35) 



(36) 



is the classical superfield = {(p) with (j) — {A, c), 
where in a slight abuse of notation we use the same sym- 
bols for the components of and (p. The flow of the 
effective action reads (cf. Eq. (fT3| ) 



2 SI6I 



(37) 



By means of ((36|) one derives from ([35]) (cf. Eq. ([14 
S^Wk fS^T 



5161 



Ilk 



(38) 



and obtains the desired flow equation (cf. Eq. (ITEl) ) 



dtVk 



-STr 



ilk (l 



.(2) 



7^^ 



(39) 



where F 



(2) 



In the present Hamiltonian 



approach to Yang-Mills theory in Coulomb gauge the ef- 
fective action Ffc[0] defined by Eqs. ^ and ^ 
is exclusively determined by the vacuum wave functional 
'ip[A\ and the Faddcev- Popov determinant. Importantly, 
the ERG approach docs not require the knowledge of the 
full vacuum wave functional. It is sufficient to know the 
wave functional in the asymptotic region A: — >■ oo, where 
perturbation theory applies. The full quantum effective 
action Pfc-j-o and hence the full vacuum wave functional 
is then computed by solving the flow equation, making 
suitable ansatze and truncations for P^ or its derivatives. 
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III. FLOWS FOR COULOMB GAUGE YANG 
MILLS THEORY 

A. Uniqueness of IR scaling 

The Hamiltonian flow equation set up for Coulomb 
gauge Yang-Mills theory in the last section IIIBI very 
much resembles the one in Landau gauge, but in one di- 
mension less. It has already been speculated that there is 
a close connection between these two formulations, j23j . 
Here we shall employ the similarities in order to derive 
unique scaling laws for the infrared behavior of Coulomb 
gauge Yang-Mills theory. 

It has been shown in [23, [1^ that Landau gauge Yang- 
Mills theory admits a unique infrared scaling solution 
[26| . Moreover, this solution implies ghost dominance. 
We emphasize that uniqueness refers to the unique scal- 
ing relations if scaling is present. Indeed, Landau gauge 
Yang-Mills theory also admits a solution without such a 
scaling behavior, the decoupling solution. More details 
can be found e.g. in The scaling and decoupling 
solutions also exist for the DSE obtained in the Hamilto- 
nian formulation of Yang- Mills theory in Coulomb gauge 
and were baptized "critical" and "subcritical" solutions 
[l9| . Furthermore, lattice calculations [l^l show that the 
scaling or critical solution is realized in Coulomb gauge. 

The proof in [23] was based on the comparison of the 
full hierarchies of DSE and FRG equations for Green's 
functions. Apart from this it only relied on the details of 
the coupling between ghosts and gluons and the canon- 
ical scaling properties of the gluonic self-coupling. This 
has been made transparent in [23 |. The proof as formu- 
lated there can be directly transferred to Coulomb gauge, 
the only missing piece is provided by the flow equation 
derived in the present work. 

With the Hamiltonian Coulomb gauge DSEs and the 
FRGs we derive the same set of constraint equations for 
the scaling coefficients as in [IJ]. There are additional 
terms coming from the higher classical gluonic n-point 
vertices which all can be proven to be sub-leading. This 
relates to the fact that the canonical scaling of classical 
gluonic vertices is less divergent than that of the dressed 
vertices. In summary we conclude that also Coulomb 
gauge in its Hamiltonian formulation admits a unique 
scaling solution with the same scaling laws that are sat- 
isfied in Landau gauge in c? = 3. The scaling relations 
relevant for the present work are that for the propagators. 



{A{p)A{-p)) oc , {c{p)c{-p)) « p2(i+K.) 



(40) 



and for the ghost gluon vertex at the symmetric point. 



oc p 



S^S 



(41) 



The scaling solution entails the non-renormalization of 
the ghost-gluon vertex, KccA = 0, and a scaling relation 



for the scaling of the ghost and gluon propagator, sum- 
marized as 



f^ccA = and KA 



2kc , with Kj^ < 



4 ' 

(42) 

Eq. (|42p implies ghost-dominance in the sense that di- 
agrams with ghost lines dominate in the infrared over 
diagrams with gluonic lines, see also [l^l- The scaling 
coefficients a,/3 used in Coulomb gauge are defined via 



(A(p)A(-p)) cxp", (c(p)c(-p)) cx 



1 



p 



2+13 ' 



(43) 



see Eq. (|58p below. The coefficients a and /3 relate to 
the k's via a = —2 — 2k.a and P = 2kc. Hence we find 
unique scaling laws in Coulomb gauge with the scaling 
relation 



2/3-1 



(44) 



The sum rule (|44p has been found in DSE analyses 
in Coulomb gauge before [l7j . here we have proven its 
uniqueness. 



B. Propagator flows 

The propagator flows are obtained from the flow equa- 
tion for the effective action (p9)) by differentiating twice 
w.r.t. the flelds. These derivations are detailed in Ap- 
pendix |b1 their outcome is represented diagrammatically 
in Figs, m El All propagators and vertices, denoted by 
black and white circles respectively, are fully dressed k- 
dependent correlation functions. This has to be com- 
pared with DSE equations for the propagators where all 
diagrams contain one bare vertex. 

For a solution of these equations we have to approx- 
imate the full effective action. In the present work we 
approximate it by the classical action and fully momen- 
tum dependent (inverse) propagators. We are specifi- 
cally interested in the infrared where we assume scaling. 
Uniqueness of scaling as proven in the last section IIII Al 
then implies ghost dominance. Consequently, we drop 
the gluonic vertices and only keep the ghost vertices. The 
resulting flow equations for the gluon and ghost propa- 
gators are shown in Figs. 21 |5l 

We pause here for a moment to discuss the meaning of 
and the justification for this truncation in detail. The 
generating functional (j^^ is a functional integral de- 
fined by the full vacuum wave functional i)[A\ which is, 
however, unknown. In Ref. [27j . the vacuum functional 
has been determined explicitly to one-loop order through 
a perturbativc solution of the Schrodinger equation for 
the Christ-Lee Hamiltonian. It was found that non-local 
terms in the couplings of two, three and four gluon fields 
in the wave functional give contributions to the static 
gluon propagator that are relevant to its ultraviolet be- 
havior, in particular its anomalous dimension. In higher 
loop orders, non-local terms in the couplings of more than 




FIG. 2: Flow equation of the gluon propagator, Eq. HB9|I . Here and in the following, the spiral and dotted lines with the black 
circles denote the regularized gluon and ghost propagators at cutoff momentum k, respectively. White circles stand for proper 
vertices at cutoff k, a regulator insertion 7?^ is represented by a square with a cross. 




FIG. 3: Flow equation of the ghost propagator, Eq. (|B8P 



four gluon fields will also become relevant to the gluon 
two-point function. By non-local we refer to coefficient 
functions that become singular for exceptional momenta. 
By neglecting these terms in the truncation considered 
in the present paper, the ultraviolet behavior of the two- 
point function will not be accurately reproduced, i.e., 
the power of the logarithmic correction in this kinematic 
regime will be incorrect. 

On the other hand, such non-local terms and three- 
and four-gluon couplings are not necessarily relevant to 
the infrared behavior, which is our main concern here: 
it has been argued [28| that in the infrared the static 
ghost propagator is strongly enhanced relative to its tree 
level behavior, while the gluon propagator is suppressed 
or even vanishing. This is precisely what happens for 
the unique scaling solution as discussed in the previ- 
ous section: the infrared behavior is dominated by those 
diagrams with the largest number of ghost propagators 
("ghost dominance"). Indeed, the arguments about the 
kinematic singularities in can also be directly trans- 
ferred to Coulomb gauge. They entail that for the scal- 
ing solution neither the non-local terms described before 
nor the couplings of three or more gluons in general con- 
tribute to those diagrams that dominate the infrared be- 
havior. The same can be inferred from the diagrammatic 
analysis of Ref. [l^l when extended to higher perturbative 
orders, if one takes into account that there are no non- 
local terms or higher couplings including ghost fields in 
the "action" ((26)) . 

The approximation of keeping a bare or tree-level 
ghost-gluon vertex is based on the "non-renormalization 
theorem" for this vertex [2^. Although this theorem 
was originally formulated for QCD in Landau gauge and 
the space-time correlation functions, the argument car- 
ries over without change to the present situation. Indeed 



we have shown in the last section that it follows for the 
unique scaling solution. It has been confirmed on the 
non-perturbative level for the Landau gauge case in lat- 
tice studies d^]. As for the Coulomb gauge, a pertur- 
bative evaluation of the vertex (to one-loop level) at the 
symmetric point shows that the quantum corrections are 
finite and independent of the scale [13, [3l| . 

In summary, we can drop the gluonic vertices in the 
infrared without spoiling the quantitative nature of our 
approximation. We emphasize that for large momenta 
this is evidently not true. Finally, we also drop the tad- 
pole diagrams in the fiow equations for the static prop- 
agators. The four-point couplings appearing in these di- 
agrams are not contained in the integrand of the gen- 
erating functional, but can build up in the course of the 
rcnormalization group fiow. We simply assume that their 
contribution is negligible in the infrared, at least for the 
qualitative behavior of the two-point correlation func- 
tions, which leads us to the final form of the fiow equa- 
tions represented in Figs. IH [S] 



IV. EFFECTIVE ACTION AND FLOWS 

A. Parametrization of the effective action 

The arguments above single out a specific approxima- 
tion of the effective action. First of all, it relies on an 
expansion of the effective action in powers of the field, 

(45) 

where the (j)^ stand for either the gluon fields {(j> = A) or 
the ghost fields {(j> = c,c). In the following we take into 
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kdk ^m^^^m^ ^ = - m5* ^ - '^(m 

FIG. 4: Truncated flow equation of the gluon propagator. Here and in the following, the bare vertices at /c = A are symbolized 
by small dots. 



kdk 




+ 




FIG. 5: Truncated flow equation of the ghost propagator. 



account the bare ghost-gluon vertex and the full momen- 
tum dependent propagators. In an upgrade of the ap- 
proximation we will also take into account tadpole terms 
related to ghost and ghost-gluon vertices. This provides 
a first error estimate for the approximation scheme set 
up here. All other vertices are set to zero. Therefore, 
in the minimal order of the approximation the only non- 
vanishing n-point functions are the ghost and gluon (in- 
verse) propagators and the ghost-gluon vertex. 

We parametrize the inverse gluon propagator as fol- 
lows, 



with 



(46) 

The diagonality in color space is due to global color sym- 
metry. The transverse projector comes with the choice 
of Coulomb gauge where the gauge fields are transverse, 
see (j2ip , and momentum conservation arises from spatial 
translational invariance of the theory. The only quan- 
tity left to be determined by the flow equation is LUk{p), 
which depends only on the absolute value of the exter- 
nal momentum due to rotational invariance of the theory, 
and on the cutoff momentum k. The factor of 2 is mere 
convention. In the flow we need the gluon propagator 
GA,k tijS"'^ with the scalar function 



GaAp) = 



1 



2uJkip) + RA.kij)) 



(47) 



The full gluon propagator at vanishing cut-off is given 
by Ga{p) = l/2a;(p), with w(p) = a;o(p). The ghost 
two-point function can be expressed as 



^'^^'"'5:T^(27r)^'^^(p + q) , (48) 
dk[p) 



Jc''(p)<5c''(q) 



where dk{p) is the ghost form factor, which is the quan- 
tity to be calculated. The ghost propagator Gck 5°''^ com- 
prises the scalar function 



Ge,fc(p) = -G,.k{p) 
9 



(49) 



GcAp) = 



P^/dk{p)+RcMip) 



(50) 



We have included an explicit constant factor of 1 /g in the 
definition of the ghost form factor for ease of comparison 
with the Dyson-Schwinger equations of the variational 
approach in subsection IIVCI The full ghost propagator 
at vanishing cut-off is Gc(p) — d{p)/gp^, where d{p) = 
do{p). The last quantity to specify is the ghost-gluon 
vertex. We have argued in the previous section that it is 
well approximated by its bare part. 



-igr'''pi,jUjiP3){27rfS^iPi +P2 + P3) , 



(51) 



where we have used the fact that the gauge field is trans- 
verse in Coulomb gauge. With these conventions, in par- 
ticular Eq. (|48|) , a suitable choice of the renormalization 
group invariant coupling gjt is 



gnip) = d{p){pMp)y/' , 



(52) 



see [32 



B. Approximation without tadpoles 

Now we are in the position to solve the flow equations 
for the propagators within the minimal truncation in- 
troduced above: the only non-vanishing vertex function 
is the bare ghost-gluon vertex ([ST|) . In particular this 
eliminates the tadpole diagrams. Inserting this approx- 
imation into the flow equations for the gluon and ghost 
propagators shown in Figs. 21 Owe arrive at 
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(53) 



cPq r 



(27r)3 



GA,kRA,kGA,k (q) Gc.k{\p + q[ 



2 

(G,.fci?,,fcGe,fc) (g) GaAIp + q|) (i - (p • q)') 

V y (p + q)^-i 



(54) 



The computation is detailed in Appendix C. Eqs. ()53p 
and (j54p are two coupled functional ordinary differential 
equations for and dk , which can be solved numerically. 
Due to our definition p8)) of the ghost form factor dk {p) 
the bare coupling constant g has formally disappeared 
from the propagator flow equations. Let us stress again 
that a physically meaningful rcnormalization group in- 
variant coupling is defined by Eq. (|52p . To incorporate 
the appropriate initial conditions it is convenient to cast 
the differential flow equations into an integral form, 

LOk{p)~u^(j>) = j^—j ItAdk'W,v) , (55) 

-v y (^4[^fcsd.'](€,p). 

(56) 

and I'^ stand for the integrands of the loop integrals 
on the r.h.s. of Eqs. (I53|) and (|54|) . 

The initial conditions d/^{p),uj/^(j)) for the flow can be 
determined by perturbation theory [l^, [HI . Due to the 
kinematic structure of the ghost-gluon vertex no mass 
term for the ghost is produced. Moreover, contributions 



with higher powers of momenta than the ones in lowest- 
order perturbation theory, d^^\p) = 1 and uj'^^'>{p) = p, 
are suppressed by the corresponding powers of A. Note, 
however, that there is additional logarithmic scaling. In 
the case of d/^(p), this implies that for a large initial cut- 
off scale fc = A we only have to fix the constant d^ip) = 
dA- 

For the gluon the introduction of the regulator term 
enforces a mass- like term, i.e., a p- independent contribu- 
tion to ujk (p) , due to the modified Slavnov- Taylor identi- 
ties (mSTI), see e.g. 0,11, [13, [Hill. We rush to add that 
the direct use of the mSTI in the infrared is not advisable 
as the mSTI only fix the longitudinal mass which does not 
relate to the transverse one in the infrared [l|, d, [11] . We 
conclude that for large initial cut-off scales fc = A the in- 
verse gluon propagator has two relevant parameters, the 
mass-like parameter a and the coefficient of the classical 
term. The latter is put to one and wc have wa(p) = P+a. 

It is evident from the form of the flow equations in 
(|55|) , ([55)1 that the solution will not show infrared scaling 
unless the parameter d~^^{p) = is fine-tuned (at least 
for /? = 2kc > 0). Such fine-tuning of relevant parameters 
is a well-known initial condition problem for RG flows. 
Indeed, it also occurs in Landau gauge Yang-Mills the- 
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k 



FIG. 7: The flow of the fuU ghost form factor, dk(p), is shown. 
The gradual formation of the IR power law when lowering the 
cutoff scale k is explicitly seen. 

ory, see e.g. [H, IH, [s^ , where it is directly related to the 
resolution of the Gribov problem, see [H, 13 • In [H it has 
also been shown that there is a full one parameter family 
of solutions compatible with non-perturbative renormal- 
ization where only the endpoint shows a scaling behavior 
whereas the other solutions show a decoupling behavior: 
a gluon with a mass-like propagator and a ghost which is 
at most logarithmically enhanced. Such a scenario very 
likely also applies to Coulomb gauge and hence deserves 
further investigation. 

In the present case we have numerically solved the fine- 
tuning condition for c?a with the constraint of infrared 
scaling for the ghost dressing function. The parameter a 
in the initial condition uja(j)) is fixed by demanding that 
oj{p) reduce to the perturbative form uj{p) (x p for "large" 
momenta p close to but below A. The regulator used in 
the numerical solution is 

rfc(p) =cxp(^^-|j^ . (57) 

Our numerical procedure is detailed in Appendix |D] The 
results for the inverse gluon propagator uik{p) and the 
ghost dressing function dk{p) arc shown in Fig. [6] for dif- 
ferent values of the minimal cutoff kmin down to which 
the flow integration has been carried out. It is seen that 
the power law behavior in both cases extends towards the 
IR as the cutoff kmin is lowered, although we have imple- 
mented a scaling behavior - not the "horizon condition" 
d-o^iP = 0) = - only for the ghost dressing. Thus, we 
arrive at the nontrivial result that a solution for the flow 
equations can be found that obeys a power law behavior 
(HSl) for both the gluon and the ghost propagator. For 
the sake of illustration, we also display the full flow of 
the ghost dressing function, dk{p), in Fig. [T] 

The IR power laws are extracted from the numerical 
solution shown in Fig. [HI The scaling coefficients a, f3 
defined in arc computed as 

w(p^O)-p~", d{p ^ 0) r--. p-P . (58) 



and with the numerical solution we get 

a = 0.28, P = 0.64 , (59) 

and hence a and (3 satisfy the sum rule already 
found for the Coulomb gauge DSE in (T^]. Note how- 
ever that the scaling coefficients obtained in the present 
truncation differ from the ones obtained in the DSE. In 
Fig. [5] the solutions of the FRG flow equation for uj{p) 
and d{p) arc compared to the results obtained from an 
optimized calculation in sect . IIV CI which in turn are 
precisely the results found in jlq bv a variational calcu- 
lation, see Fig. IHl This variational calculation gave rise 
to the DSEs which will be found in sect. IIV CI as approxi- 
mation to the full flow equation. While the curves in Fig. 
[5] match in the UV, the results of the FRG in the present 
minimal truncation arc less infrared enhanced than the 
ones of the DSE. Note that the scaling coefficients are 
expected to depend on the chosen regulator. It has been 
already proven in [s^ for Landau gauge Yang-Mills the- 
ory that the scaling coefficients of FRG and DSE agree 
for optimized regulators if a bare ghost-gluon vertex is 
used. Optimized regulators are those that maximize the 
physics content of a given truncation scheme at a given 
order, for details of the optimization theory for the FRG 
we refer the reader to la. [371. The details of the ar- 
guments put forward in jSa] directly carry over to the 
Coulomb gauge system. 

C. Optimization 

In this chapter we shall use the optimization arguments 
hinted at above in order to optimize the physics content 
of our present truncation. Similar, but more refined ar- 
guments have been used in Landau gauge for arriving 
at FRG results for the propagators that quantitatively 
agree with the lattice results [l|, [35| , the details will be 
published in |35| . Here we first follow the arguments 
put forward in [36|. To that end consider the following 
approximation: under the loop integrals we replace the 
propagators at the running momentum scale k,ujk and 
dk, by the propagators of the full theory, i.e. the ones at 
zero scale fc = 0: 

dkip) dk=G{p) , iOk{p) ^^k=o{p) ■ (60) 

Note that this does only imply that the difference be- 
tween the propagators at fc = and the regularized ones 
at k drops out in the integrals. Indeed one can explic- 
itly construct regulators for which this holds true in the 
asymptotic IR region, see [s^. Note that due to the 
strong infrared suppression introduced with the regulator 
choice (j57p the approximation (|60p is quantitatively reli- 
able inside the loop integrals except for a small range of 
momenta p around the scale k. The approximation (j60p 
allows us to analytically integrate the flow equations (|53p , 
(IMl) over k. The only fc-dependence left is the explicit 
one on the regulator as the vertices are not fc-dependcnt 
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from the outset. Hence the flow can be rewritten as a to- and we arrive at 
tal t-derivative of the loop integrals with full propagators 



{ujQ - uja){p) 



-N, 



1 



(i-(p-qn 



(27r)3 2^o(g) + RaAq) (p + q)2do ^(|p + q|) + Rc,ki\P + ^1 
Nr. r d'q (l-(p-q)') 



1 



4 J (27r)3 g2rf-l(g)+i?,,fc(g) (p + q)2d-i(|p + q|) + i?,_,(|p + q 



fc=0 

fc=A 
fc=0 

fe=A 



(61) 



(62) 



We notice that the flow equations (|6T|) . (|62|) have ac- 
quired the form of DSEs. Given the fact that Rk=o = 
(see Eq. (|3]) these equations coincide with the DSEs 
obtained in [l^ (to be precise, the DSE for uj(p) in 
[isj contains additional contributions which are, how- 
ever, subleading in the infrared), with a different UV- 
regularization realized here through the (fc = A)-terms. 
Moreover, the optimization arguments in 0, [36| imply 
that the flows (|5T|) . provide the best approximation 
to the full theory for the IR asymptotics. In [s^ it is 
shown that in the given truncation this argument extends 
to the full momentum regime: by adding the tadpole di- 
agrams related to ghost-only and ghost-gluon vertices to 
the flow equations displayed in Figs. HI [5] and using the 
DSE for the tadpole vertices in the flows one can show 
that this leads to the integrated flow (|6ip . If we 

also add the gluonic diagrams, this argument gets more 
involved. 

It remains to adjust the initial conditions. We could 
proceed in the same way as for the the numerical solution 
in the last subsection to implement the condition of in- 
frared scaling for the ghost dressing function. However, 



it is much simpler to use as an input the information 
from this numerical solution that /3 = 2k > and thus 
d^'^{p = Q) =Q (the horizon condition), so we can write 



[mt{k, p, q) - int(fc, p = 0, q)] 



fe=0 



fc=A 

(63) 

where int(fc, p, q) denotes the integrand in (|6T|) . More 
details of the numerical procedure can be found in Ap- 
pendix |D] 

The results of the iteration are shown in Fig. [9l A 
power law as in Eq. (|58p emerges in the infrared region 
for both the gluon energy uj{p) and the ghost dressing 
function d{p) with the IR exponents 

a = 0.60, P = 0.80 , (64) 

which is precisely one of the two possible IR solutions 



found analytically for the DSE in [l7[. Furthermore, 
it corresponds to the solution found in the variational 
approach [l^. It must be noted that the second pos- 
sible solution in the analytical approach of Ref. [l3l has 
also been confirmed in a numerical variational calculation 
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gluon correlator CO (FRG) 

— ghost dressing d (FRG) 
• ghost dressing d (DSE) 
■ gluon correlator CO (DSE) 




FIG. 9: Gluon propagator lj and ghost dressing function d 
from the optimized flow equation in comparison with the re- 
sults of the DSEs obtained from the variational ansatz in [3| • 
The results lie on top of each other as expected. 
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FIG. 10: The running coupling constant a = gji/4iT, Eq. ([52}, 
calculated from the optimized flow and from the flow without 
tadpoles. 



|18| . Indeed, it is also present in our optimized approxi- 
mation which has the identical infrared properties as the 
DSE. However, it is not clear to us at present whether 
after the inclusion of the gluonic diagrams this solution 
persists as an infrared stable one. Note in this context 
that it requires additional fine-tuning and hence may be 
unstable. 

Finally, in Fig. I10| wc show the running coupling con- 
stant a = gj^/An, see Eq. ([52|) . calculated by using the 
propagators from the optimized flow equation as well as 
from the flow equation without the tadpoles. The plateau 
in the IR is due to the sum rule (|44l) which is fulfilled by 
the propagators resulting from both approximations of 
the fiow equations. 



V. CONCLUSIONS 

We have presented a new approach to the non- 
perturbative calculation of static propagators, or corre- 
lation functions at equal times, from the Hamiltonian 
formulation of a quantum field theory. In the generat- 
ing functional of the correlation functions at equal times, 
(minus twice) the logarithm of the modulus of the vac- 
uum wave functional comes to play the role of the Eu- 
clidean action in the usual covariant theory We have 
then adapted the functional rcnormalization group to this 
Hamiltonian formulation. 

This new approach has subsequently been applied to 
Yang-Mills theory in the Coulomb gauge. The derivation 
of the corresponding flow equations for the propagators 
has been presented in detail. In order to arrive at a closed 
system of equations, we have replaced the dressed ghost- 
gluon vertex with the corresponding bare one (or a vertex 
with constant dressing, a good approximation according 
to perturbative arguments and lattice calculations), and 
we have neglected all tadpole diagrams and vertices with 
three or more gluon lines. We have also presented an 
approximation that allows for an analytical integration 
of the flow equations and have argued that it actually 
corresponds to an optimized choice of the regulator func- 
tions. We have discussed in detail the choice of the initial 
conditions which serve to implement the normalization 
conditions corresponding to an infrared scaling solution 
for the propagators. 

The result of the numerical solution of the flow equa- 
tions has been compared to the solution of a system 
of Dyson-Schwingcr equations derived from a variational 
principle for the vacuum wave functional. The solution 
of the flow equations agrees with one of the two possi- 
ble scaling solutions in the latter approach. In the ap- 
proximation without tadpoles we have obtained slightly 
different values for the infrared anomalous dimensions 
than from the Dyson-Schwinger equations. Indeed, in 
this truncation the anomalous dimensions mildly depend 
on the regulator functions and only agree with the Dyson- 
Schwinger values for optimized regulators, as has been 
argued in [36j . In the optimized approximation the in- 
frared anomalous dimensions are regulator-independent. 
This supports the reliability of the optimized approxima- 
tion. 

Wc have not presented the second scaling solution 
which poses an additional fine-tuning problem. Its res- 
olution also allows to study the interesting question of 
infrared stability of this solution, and will be discussed 
elsewhere. 

Note also that the formalism put forward in the present 
work allows to directly access the confining properties of 
the propagators. This can be done by using the Wilson 
loop potential evaluated in [38[ . Moreover one can study 
full dynamical QCD in Coulomb gauge along the lines of 
[39| . A first step in this direction is to implement the 
missing gluonic diagrams for a comparison with lattice 
data in the full momentum regime. These investigations 
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will certainly shed more light on the pressing unresolved 
questions of low energy QCD. 
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Appendix A: Details of the derivation of the flow of 
the effective action 



which turns into the usual one upon taking fc — )■ 0, be- 
cause i?A,fc=o = Rc,k=o = 0. In Eq. (|A4p . the sources are 
functionals of the fields (whose notation is suppressed), 
which are the expectation values of the corresponding 
field operators. The relations between sources and fields 
are given by 



^n-p)=(27r) 
c"(-p)=(27r 



^ 5Wk[Jk,<Jk,<Jk\ 
(5ct°(p) 



(27r) 



5(7" (p) 



(A5) 



In this appendix we fill in some of the details of the 
derivation of the flow equation ((39|) . From the regularized 
generating functional (^5)) and Eq. (P7)) . we derive the 
flow equation 



j^-Rc.k--h)Zk[Ji<yi<y] , 
(Ai) 

where the difference in signs of the gluon and ghost regu- 
lator term compared to Eq. p7)) is due to the Grassmann 
property of the ghost fields and sources. 

The definition of the Schwinger functional Wk gener- 
ating connected Green functions reads 



Wk [J, cr, cr] := In Zk [J, cr, cr] , 
and therefore its fiow is 



(A2) 



dtWk 



l5Wk ■ 5Wk 

tiA.k 



2 5.J 

5Wk 



5J 



Sa 



■R. 



c,k 



SWk 
5a 



TtR 



■c.k 



SaSa 



(A3) 



The effective action is defined via a modified Legendre 
transformation. 



rfc[A, c, c] -.^ -Wk[Jk,o-k,ak] + Jk-A + ak-c + c- (Jk 
- ■ RA,k ■ A-c- i?c,fe ■ c , 

(A4) 



or shorthand 



A" 



5Wk 



5a ' 



5Wk 
' 5a 



(A6) 



With these definitions and with Eq. (|A3|) , the flow of the 
effective action is written as 



't-L k 



{dtWk)[Jk,Ok,(^k] 

^ J SWk „ SWk ^ _ SWk 
• —J Otak ■ Otak ■ -jT- 

A ■ dtJk + dtak ■ c - dtak ■ c 
^A ■ ilA^k ■ A ~ c ■ Rc,k ■ c 

5^Wk 



1^ p ^'^fc 



Tr Rc.k 



SaSa 



(A7) 



When expressed in the superfield notation introduced 
from Eq. ([3T|) on, this equation turns into Eq. (f37| and 
finally into Eq. ([M)) . In components of the superfield, 
Eq. §9^ reads 
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dtTk 



■Tr 



T 

'c,k 



5A5A 



RA.k 



Sc6A 

ScSA 
I 



SASc 



dcSc 



R. 



■c.k 



ScSc 



SASc 
ScSc 



ScSc 



(A8) 



From this equation wc derive the flow equations for the 
ghost and gluon fields by taking functional derivatives 
w.r.t. these fields. 



Appendix B: Details of the derivation of the 
propagator flows 

In this section we derive the flow equations for the 
propagators from the flow equation for the effective ac- 
tion, Eq. (|A8[) . A useful relation for the following con- 
cerns the commutation of fermionic derivatives past su- 
permatrices. Consider a matrix with a block structure of 
commuting (c) and anticommuting (a) quantities and an 



anticommuting 77: 





a 








c 








c 










a 


a\ 


]•-'{: 


c 


c \ M r] 




c 





(Bl) 

This matrix structure is shared by S'^T^/ S4>5(j) as well as 
{5'^Tk/5(t)S(l) + 7^fc)~^ the latter because of Eq. (pS)) . 

In the following, i and j are condensed external indices. 
They stand for color indices, momenta and, in the case of 
gluon fields, also vector indices at the same time. They 
are, however, not part of the matrix notation. From 



J 



it follows that 



= 



Sci \ \ S(j)S4) 



SciScjyScf) \S(j)S(j) 



M 



Uk M- 



S fS^Tk 



-R-k 



Sci \S(l)S<j) 



S^'Tk fS^Tk 



\S(f)S(l) ^) SciS(f)S(f) \S(j>S4> 



7^^ 



(B2) 



(B3) 



Sci \64>S(f) 

For bosonic derivatives the same formula holds without the M's. Using this, we can derive the ghost propagator flow 
equation from the flow of the effective action, Eq. (jASP : 



S'^Tk 1 S'- 



ScjSci 2 ScjSci 



STr 



Ilk 



1_J_ 

4 STr 



STr 



-7^fc M 



Ilk M 



S^Vk fS^Vk 



SciS(l)S(j) \S(j)S(j) 



TZk M- 



ScjScjjScj) \S(f>5(j) 



-1 

7^fc 1 M 



SciS(j)S4) \S(j)S(f> 



7^^ 



(B4) 



-STr 



S^T 



TZk M- 



SciS(j)S(f> \S(j)5(j) 



Tlk M 



ScjS(j)S(j) \S(f>S(j) 



7^fe 



- -STr 
2 



s'n 
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Setting the fields to zero, A 
remain. With the definition 



0, only the block matrices with the same number of ghosts and antighosts 



6A6A 



RA,k 



ScSc 



the first of the three terms in Eq. (jB4p becomes 










\ 


1 





— Rc,k 







I 





~^c,kj 



Qk 









Scj6A6c 



\ ScjScSA 



\ ScSc 



Rc.k 







SciScSA 




SciSASc 




Tr RA,k 



6A6A 



1 



■ Tr i?c,fc 



ScSc 



RA,k 



^c.k 



ScjSASc 



( 



ScSc 



S'T, 



ScjScSA \SASA 



Rc 



Ra.. 



\ 







Qk 



J 



ScjScSA \SASA 
S^Tu fS^Vk 



SciSASc V ScSc 



Ra,^ 



RJl 



The other two terms are treated alike. The ghost flow equation then reads (using f, = Rc.k) 



ScjSci 



=Tr RA.k 



( <s^r, 



\SASA 



TtR 



■c,k 



ScSc 



Ra., 



+ Rc.k 



S'T, 



ScjSASc 



ScSc 



-1 



SASA 



Ra., 



Sc^ScSA \SASA 

-1 c- 



+ Rc.k 



Ra., 



s'r, 



fn_ 

SciScSA I SASA 



Ra.i 



S^Tk 
ScjSASc 



ScSc 



+ Rc.k 



fS^Tk 



ScjSciSASA \SASA 



+ RA.k 



TtRc 



ScSc 



Rc 



s^r. 



ScjSciScSc 



ScSc 



R. 



In much the same way, the ghion flow equation is deduced from Eq. (jASp . and the result is 



SAjSA, 



=Tr Ra., 



TtR, 



TtR, 



fs^r, 



\SASA 



Ra.i 



S'Tk f S^T, 



~.k 



~.k 



■ Tr RA.^ 



ScSc 

S'Tk 
ScSc 



Rc.k 



SAjSASA \SASA 



S^T, 



SAjScSc \ ScSc 



RA.k 



Rc 



S^T, 



SA.SASA \SASA 



s'r. 



SAiScSc \ ScSc 



S'Tk 



R. 



S^T, 



\SASA 



c.k 



Ra,, 



SAiScSc 
1 



S'T, 



+ R. 



ScSc 



c.k 



SAjScSc 



ScSc 



Ra., 



+ Rc.k 



+ Rc.k 



SAjSA,SASA \SASA 



R 



A.k 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



- Tr R, 



C.k 



ScSc 



R. 



c.k 



k 



SAjSAiScSc 



ScSc 



R, 



c,k 



These equations are represented diagrammatically in Figs. [21 El 
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Appendix C: Explicit form of the flow equations 



111 this appendix we will use the parametrization of sect. IIV Al in order to bring the truncated flow equations 
represented in Figs. 21 [S] into their final form. 

As the ghost two-point function and the ghost regulator are both diagonal in color space and momentum space, 
they can easily be inverted, yielding (see Eqs. ([29]), (gS]), dH]), dSO]) ) 



6^ 
ScSc 



+ Rc.k 



ah 



S""' [gpVdkip) + RcAp)] ' i^^fs'ip + q) 



9 



(CI) 



A similar formula holds for the gluon two-point function and the gluon regulator, which are invertible in the transverse 
subspace (see Eqs. ([211), (05]), (|T7)) ): 



ah 



= (5"%- (p)GA,fe(p) (2^)^(5^ (p + q) 



Therefore, the truncated gluon flow equation. Fig. 21 see also Eq. (|B9|) . reduces to 
2<5%(p)c:;fc(p)(2^)353(p + q) = 

<5'=^4,fc(pi)(27r)3<53(pi - p2) 5'''G,M{2nf5^{Y>2 - Pa) 

X rfUM{ip4,n.)i27Tf6^p + P3 - P4) 5fSG,M{p4){27TfSHp4 - Ps) 

X fs''H,M){tpej){2nf5^{q + Ps - Pe) S'^-^G^Mi^T^fS^iPe - Pi) + (« ^ j) 
and by carrying out index contractions and integrations we obtain the flow equation for ujk'. 



(27r)3 



Gc,kRcMGc,k) {q)GcM\p + q|) ^'(i - (p ■ q)') 



In much the same way we treat the truncated ghost flow equation. Fig. [51 sec also Eq. (|B8p . to get 
-rtd^\p)(27r)3<53(p + q) = 



rf'[pi...6] 
(27r)i8 



<5^%(pi)i?A,fe(pi)(2^)3<5^(pi - P2) 5^H,hiP2)GA,kiP2)i'2nrS'(.P2 - Pa) 
X /^^'^ti^(p3)(-iP4,/)(27r)3(53(p + p3 _ p^) SfsG,,k{P4){27TfS^{P4 - P5) 
X /"'^tm„(P6)(*gm)(27r)3,53(q + p5 -pe) 6''^tm{pe)GA,k{Pe)i'2^f^HP6 - Pi) 

S'^'^R^MPlK^T^fS^Pl - P2) S^''G,.k{p2){27Tf5\p2 - P3) 



X f'\m{Pi){iqrn){2TTfS''{q+p3 - P4) <5'=-^ty (p4)GA.fe (^4) (2^)^(53 (p4 - Ps 

X rf%,j{p,)Hpe,h){27rfS^{p + P5 - Pe) <5^''Ge,fc(p6)(27r)3<53(p6 _ pi) , 
which, after receiving the same treatment as the gluon equation above, finally turns into, 

d^q 



dtd,\p)^N,p' 



J (^ {GA,kRA,kGA,k) iq)Gc,ki\p + q|)(i - (p ■ q)') 



{Gc,kRc,kGc,kj iq)GA,ki 



p + q|) 9 



1 - (P ■ q)^ 

(p + q)2 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 
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Appendix D: Numerical Computation 

Below we give some details of the numerical method 
used for the solution of the flow equations. In order to 
solve the flow equations numerically, the functions u}(j>) 
and d{p) of section llV Cl as well as uJk{p) and dk{p) of sec- 
tion IIVBI are represented by means of Chebyshev poly- 
nomials on suitably sized momentum ranges, spanning 
up to nine orders of magnitude. We use a logarithmic 
momentum scale to sample the behavior of the functions 
equally well on all orders of magnitude; the function val- 
ues are also represented logarithmically. In the case of 
the flow functions u>k{p) and dk{p), which possess two 
momentum arguments, we first calculate the coefficients 
along the p-direction at constant values of k which are 
the Chebyshev nodes of the fc-range, then determine the 
coefficients of the resulting functions in fc-direction for 
each p-coefficient. About 130 Chebyshev nodes in each 
direction have been used. The momentum integrals, for 
the loop as well as for the flow (in section |IVB[) . have 
been calculated using the Gauss-Legcndre method with 
about 70 nodes. As for the Chebyshev nodes, the Gauss- 
Legendre nodes for the momentum integrals have been 
calculated for a logarithmic momentum scale in order to 
appropriately sample the IR behavior. The loop integrals 
have been confined to a range of one order of magnitude 
around the cutoff momentum k outside of which there 
is virtually no contribution to the integral, owing to the 
regulator functions employed, sec Eq. (j57p . 

In the loop integrals, the flow functions need to be eval- 
uated also outside the Chebyshev representation range 
of the external momentum p. But the IR extrapola- 
tion of the flow functions requires no extra assumptions 
about their IR behavior: as long as fc„im never reaches 
the lower boundary of the p-range, the functions can be 
simply chosen to be constant beyond this boundary, as 
clearly seen in Fig. [T] For their continuation in the UV, 
a power law has been fitted to their behavior in the UV 
region of the representation range, which is U!k{p) ^ p 
and dk{p) const. As for the fc-range, the functions are 
never evaluated outside the Chebyshev range anyway, so 
no extrapolation is necessary. 

The sets of equations as well as ((5^ 

have been evaluated iterativcly: starting from constant 
functions lo and d, the r.h.s. of the equations arc cal- 
culated; from the result, the initial values of the flow, 
dA and wa, are determined following the procedures de- 



scribed in the next paragraph, giving temporary results 
^tmp ^ ^tmp^ rpp improve the convergence behavior, we 
use a relaxation method to determine the final result of 
the n-th iteration as a;^"'(p) = rw^™^(p)-|-(l— 'r)a;^"~"^^(p) 
for each Chebyshev node {k,p), likewise for dk{p)- Values 
of r e [0.1,0.5] have been used, depending on how much 
the functions change from one iteration to the next. This 
yields the new functions w^"^ and d'"' which are then 
fed back into the r.h.s. of the equations as input to the 
(n-l- l)-th step of the iteration. This iteration is repeated 
until convergence is reached. 

To determine the constant dA{p) = dA in Eq. ([55)1 . we 
demand that dfc„j„ (p) fulfill a power law for p in the IR, 
dk\ (p) -P^- To do so, we observe that for a power law 
f{p) — the expression 

P^\nf{p)^f3 (Dl) 
dp 

yields a constant value, i.e., the exponent. Let g{p) be 
the r.h.s. of Eq. ((56)) (with k = kmin)- We demand that 

|;(p|l„(,,(p) + d-))i0 (D2) 

and solve for d^^. A least squares fit of a constant func- 
tion to the expression obtained (which will in general not 
yet be constant) in the IR gives the optimal value for d^^ 
in order to achieve a power law behavior for d^^ , (p) in 
the IR region. This is done during each iteration step 
as described in the previous paragraph. In this way, we 
can impose a power law behavior on the ghost form factor 
without dictating its exponent. Concerning the Eqs. ([55)1 
(with k = kmin) and (|62[) . we fit the expression —a + bp 
to the r.h.s. in the UV and use u!a{p) = a + p to achieve 
^fe=fe™„.(p)lp^A ~ P for Eq. ([55]) and ajo(p)|p^A ^ p for 
Eq. (j62p which is an expression of asymptotic freedom. 
For Eq. ([61]) the determination of the initial conditions 
is not a numerical issue anyway, see Eq. ()63p . 

This method allows for a systematic, simultaneous de- 
termination of the solution together with originally un- 
known initial conditions to accomplish that the solution 
at fc = kmin fulfill certain properties. Solving the differ- 
ential equations in the differential form would require a 
trial-and-crror search for the initial conditions to find a 
solution with the specified properties. 
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